IMPACT  LOADING  OF  SUBMARINE  HULLS 


By 

Frederick  J.  Dzialo 

Contract  No.  0NR-N00014-68-A-0146-8 
Report  No.  THEMIS-UM-70-2 


DISTRIBUTION  OF  THIS  DOCUMENT  IS  UNLIMITED 


February  1970 


Approved  for  Release 


Reproduction  in  whole  or  in  part  is  permitted  for 
any  purpose  of  the  United  States  Government.  This 
research  was  sponsored  by  the  Office  of  Naval 
Research  under  ONR  Contract  No.  N00014-68-A-0146, 
Subcontract  No.  N00014-68-A-0146-8,  ONR  Contract 
Authority  Identification  No.  NR  200-016. 


Charles  £.  Hutchinson 
Co-Manager,  Project  THEMIS 
University  of  Massachusetts 


Ill 


AliSTUAlVT 

Following  Flugge's  exact  derivation  for  the  buckling  of  cylindrical 
shells,  the  equations  of  motion  for  dynamic  loading  of  cylindrical  shells 
subjected  to  hydrostatic  and  axial  pressure  have  been  formulated. 

The  equations  of  motion  are  applicable  for  long,  short,  or  thick 
shells,  and  are  very  useful  in  calculating  deflections  and  stresses  when 
the  impact  loads  are  applied  to  comparati vely  small  regions  of  the  shell. 
The  normal  mode  theory  was  utilized  to  provide  dynamic  solutions  for  the 
equations  of  motion. 

Solutions  are  also  provided  for  the  Timoshenko- type  theory,  and 
comparisons  are  made  between  tiie  two  theories  by  considering  and  neglect¬ 
ing  in-plane  inertia  forces. 

Comparison  of  results  is  exemplified  by  a  numerical  example  which 
considers  the  effect  of  hydrostatic  pressure  on  the  dynamic  response  of 
a  shell  simply  supported  by  a  thin  diaphragm  and  subjected  to  a  localized 
unit  radial  impulse. 
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IMPACT  LOADING  OF  SUBMARINE  HULLS 


Introduction 

The  purpose  of  this  investigation  is  to  determine  the  effect  of 
impact  leads  on  submarine  hulls.  For  ductile  materials  current  design 
methods  utilize  static  loads  for  design  with  a  performance  criterion 
that  the  hull  behave  in  a  ductile  manner  when  subjected  to  an  impact 
load  that  may  occur  during  ground  collision  or  depth  charges.  Since 
glass,  considered  as  a  possible  material  for  design  has  brittle  properties 
the  present  design  practices  must  be  re-evaluated.  Hence,  more  rigorous 
analyses  must  be  made  to  determine  the  dynamic  stress  and  deformation 
characteristics  of  glass  hulls  subjected  to  impact  loads,  and  deep  hydro¬ 
static  pressures.  The  present  investigation  will  consider  only  the 
response  of  the  shell.  The  coupled  response  of  the  shell  and  ring 
frames  will  be  considered  in  a  future  investigation. 

The  model  to  be  investigated  will  essentially  be  a  cylindrical 
shell  under  deep  hydrostatic  and  axial  pressure,  and  subjected  to  an 
impact  load. 

Following  Flugge's  exact  derivation  for  the  buckling  of  cylindrical 
shells,  the  equations  of  motion  are  formulated.  The  equations  are 
applicable  for  long,  short  or  thick  shells,  and  are  very  useful  in 
calculating  deflections  and  stresses  when  the  impact  loads  are  applied  to 
comparatively  small  regions  of  the  shell.  The  normal  mode  theory  is 
utilized  to  provide  dynamic  solutions  for  the  equations  of  motion. 
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Solutions  are  also  provided  for  the  Timoshenko-type  theory,  and 
comparisons  are  made  bctuct-u  the  two  theories  by  considering  and 
neglecting  in-plane  inertia  forces. 


Lquations  of  ilotion 


Following  Flugge's  [1]  exact  derivation  for  the  buckling  of  cylindri¬ 
cal  snells,  the  differential  equations  of  motion  for  impact  loading  of 
cylindrical  shells  under  hydrostatic  pressure  become: 


air  ♦  a..‘x  -  pa(u“  -  u'j  -  Pu"  +  a2px  =  Jia2  -~r 
a.«;  +  ad"  -  aQ  -  pa(v  +  w*)  -  Pv"  +  a2p  =  cha2  ~ r 

w  X,  )  dt' 


-ay*  -  ay"  -  aii  -  pa(u"  -  v*  +  vr)  -  Pw"  +  a V 
Z  x  f  r 


jhere 


,  .,,2  A i 

=  :na 


il*  ♦  IP 

Q.  s  _i__ — 


ir  +  ir 
q  =-i - H 
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ii  =  ^  (V*  +  w  +  .'U* )  +  7-r  (w  +  vr) 

Z  o  a 


*1  “  . 
x  a 


-  ^  (u*  +  .V*  +  M)  -  ~*r  W" 


0) 

(2) 

(3) 

(4) 

(b) 

(6) 

(7) 


(8) 


3 


)  +  |r:Li-!i  f*’  -»-) 

(9) 

'V 

!  “7  (w  +  w*‘  +  vw") 

a 

(10) 

V 

:  Tjr  (W*  +  vW“  -  U* 
a 

-  vV*) 

(11) 

"♦X 

-  jr  (1  -  »)(«’•  ♦ 

JU*  -  \  V') 

(12) 

Hx, 

=  p-(1  -  »)(«••  - 

V‘) 

(13) 

Substitution  of  equations  (4)  through  (13)  into  equations  (1)  through 
(3)  yields: 


where 


u"  +  I*  ~  v)  u“  +  ^  v'*  +  vw*  +  u*  -  w"'  +  *■••) 

-  q-j (u“  -  w')  -  q2u"  4  ^lX,bt)d  =  (14) 

Lpi  u'*  +  v  +  1 ^  v"  +  w*  +  k(  §  (1  -  v)v"  -  i^-w-) 


-  q 

vu'  +  v’  +  w  +  k 


(V  .  w.)  ■  q  Y.  +  VX»  t)a  -  Pha2  a2v 
1*V  +  Wj-q2v  +-J  g “IT  at7 


rw 

-1— “■  v  u'”  -  u"'  --3  ^  v  V"’  +  w"" 


(15) 


+  2w""  +  w"  +  2w"  +  w 


+  q^u'  -  v’  +  w) 


.  .  pr(x>  ‘)a2  .  -oha2  a*w 
q2W  D - D  W 


(16) 


k  3  ife 


ql=f 
q2  =  £ 
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Equations  (14)  through  (16)  may  be  written: 


ilr'J  +  a. 


V  '  “2 


ir  ♦  1-p  *••  +  a3W  ♦  k[lp.  - 


px(x.  t)  2  hjk2  32u 


a  = 


D  W 


u'*  +  o4(v"  +  W)  +  aiv"  -  |  (3  -  v)w- 

Pt(».  t)a2  _  ^2  s2y 

D  D  jF 

a3U*  +  a4V*  +  (2k  +  1  -  a4)W  +  kj— u'~  -  u* 
-  2-1-*  v“*  +  w“"+  2wM*‘  +  W-  +  w| 


pjx,  t)a 


where 


+  (1  -  a5)w“  - 


«,  “  i-p-  (1  ♦  3k)  -  q2 
«2  *  (1  +  k)  -  q, 


_  -pha  d;,w 
D~  at7 


<»3  =  (v  +  q^) 


,JA  =  1  -  <*1 


u5  =  1  -  q2 


Orthogonality  and  Modal  Vibrations 


For  free  vibrations,  equations  (1)  through  (3)  become 
au;  +  al4*x  -  pa(u"  -  w')  -  Pu"  =  pha2  0 
an;  +  an;,  -  aqA  -  pa(v"  +  w*)  -  Pv”  =  pha2  0 


$  X4> 


(17) 


(18) 


(19) 


(20) 

(21) 


Lquatious  (2u)  uirougii  (22)  yield  tlie  free  vibration  fmjuencies  and 
mode  shapes.  The  orthogonality  condition  is  derived  by  assuming  that  the 
displacements  u,  v,  and  w  have  the  form 


Finding  me  orthogonality  condition  involves  the  following  steps: 

(1)  the  ntn  terms  of  expressions  (23)  are  inserted  into  equations  (20) 
through  (22),  and  the  resulting  equations  are  multiplied  by  um(x),  vn(x) 
and  w  (x),  respectively,  integrated  over  the  domain,  and  added;  (2)  the 
mth  terms  of  expressions  (23)  are  inserted  into  equations  (20)  through 
(22),  and  tne  resulting  equations  are  multiplied  by  un(x),  vn(x),  and 
wn(x),  respectively,  integrated  over  tlie  domain,  and  added;  (3)  the  two 
equations  resulting  from  Step  2  are  subtracted  from  those  resulting  from 
Step  1;  they  are  integrated  by  parts,  and  use  is  made  of  equations  (4) 
through  (13)  to  obtain  the  final  orthogonality  relation.  The  orthogonality 
condition  may  be  written  as  follows: 


i  >  rL 

K  -  O  I  oh(unum  +  v  v  +  ww  )  dx 


n  m 


n  m  n  m  n  m 


Du  M  „  dv 

u  (II  +  pw  -  P  —2)  +  v  (II  -  -^2-  -  p  — £) 

n'  xm  '  m  dx  ’  V  x#i  a  dx  ' 

DW  DW 

-  w  (Q  +  P  —2)  +  II  — - 

Vwxm  dx  '  xm  dx 

DU  H  DV 

-  u  (ll  +  pw  -  P  — — )  —  v  ( N  -  -  P  — 2-) 

V  xn  1  n  dx  '  V  x«n  a  dx  ' 

DW  DW  ■>£ 

+  w  (q  +  P  --)  -  H  =0,  pi  f  n 

niwxn  dX  ’  xn  ox  L 


v 


(24) 
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where  the  natural  boundary  conditions  for  the  fixed,  simply  supported  and 
free  condition  are  given  as: 

Fixed  at  x  =  0,  ; 

u  =  v  =  w  -  =  0  (25) 

Hinge  at  x  =  0,  » 

u  =  v  =  w  =  0- 

Hx  =  C  (26) 


Simply  Supported 


at  x  =  0,  l 
w  =  0 


H  =  0 
x 


U  -  ^  -  P  =  0 
x$  a  bx 


N  ♦  pW  -  P  rj  *  0 

X  '.»X 


Free 


at  x  =  0,  •• 


v 


Xv 


jy 

a 


P  =  0 

9X 


N  +  pw  -  P  =  0 

X  oX 


(27) 


(28) 


The  differential  equations  (17)  -  (19)  may  be  solved  by  assuming 
u  =  /\e/X^a  cos(m$>) 

v  =  beu/a  sin(m<j>)  e1uit 

w  =  CeAX/^a  cos(n<j>)  e^wt 


(29) 


Inserting  equation  (29)  into  equations  (17)  -  (19)  fields  equation  (30). 


s>«^  ,L-n2^+(pl»a2i.2i«n/b)  (~g~)wA 


I 

lotgA-l: 


(-g~-)m2A-F  i"5  j  j A  j 


a '+(pha2w2nB/D)  -«^»Hc(^-)pu2  0 

a^n*k(-j-)fl).  -(2k+l-a^)n2  ♦l  C 


+(l-a5)>/ 


:  !  I 


♦k[(A2-n2)2+l]j  !  | 


•pha2u»^/u 


J 


The  characteristic  equation  is  found  by  setting  the  determinant  of 

equation  (30)  equal  to  zero.  To  determine  the  eigenvalues,  w  ,  the 

following  method  is  utilized:  A  value  of  is  guessed  and  inserted  into 

the  characteristic  equation.  The  characteristic  equation  will  yield 

eight  roots.  For  unequal  roots,  equations  (29)  may  be  written  as  follows: 

8  *.x/a  .  .  8  A.x/a  .  .. 

u  =  l  l\-a  1  (cos  m^e1a,t;) ,  v  =  7  B.e  1  (sin  ntye1  *) 
i=l  1  i=l  1 

(31) 

8  A.x/a  .  . 

w  =  l  C.e  (cos  m<f>)  e 
i=l  1 

where  for  each  A.  there  exists  a  relationship  between  the  amplitudes 
and  from  tiie  determinant  of  equation  (30). 

Equations  (31)  with  the.  necessary  boundary  conditions  will  lead  to  a 
determinant  |a- .|.  A  plot  is  then  made  of  the  determinant  j  a . - 1  versus 

1  J  I  J 

,  2  9 

</\  The  eigenvalues,  w  ,  are  those  for  which  |a.  .  |  =  0.  At  a  point,  w  , 

*  j 

when  ja^l  =  0,  the  ratio  of  the  amplitudes  ,  11.  and  can  be  calculated 
from  the  determinant  of  equation  (3u). 
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For  inpact  loads,  local  bonding  action  will  predominate,  and  trie 
principal  -;odc*  of  response  will  be  in  tnu  radial  direction.  fieglectinn 
inertia  forc-s  irr  tic  longitudinal  and  circunferer.tial  directions, 
equation  ( J'.i)  reduces  to  equation  (3/). 


'  .2 

j!V  *n  “i 

■  J  c. 

1 

(^)n. 

J  C  * 

1 

*  - 

■* 

E 

f 

H-t4"- 

-  i4m*k(-^-)m.2  | 

i 

!u‘ 

1  ' 

t 

k 

J  i  / 1  —  -  \  4,  «  3 

j  0 

; 

-(2k+l-ujm2+i  ! 

4  ‘ 

id- 

:  i 

t 

* 

+(1-^)Z  j 

I  i 

1  j 

• 

i 

f 

+fc[(*2-riz)Z+l]| 

i 

i  1 
\ 

1  i 

i  r 

i  1 

Solutions  for  Forced 

Vibrations 

(32) 


Equations  (17)  through  (19)  nay  be  solved  by  assuming 


u  =  T  u  (x,  :)  q  (t) 
n=0  n  " 


v  =  i  vn(x’  9> 
n=0 


(33) 


w  =  I  w  (x,  c)  q  ( t) 
n=C) 

Substituting  the  above  equations  into  equations  (17)  through  (19),  and 
utilizing  the  orthogonality  condition  (24)  yields  the  following: 


•2- 


r>,  ri 


1  I  [Py(x,i,  )u  +P  (x,^,a)v  +P  (x,ij),Ov^  ][sin  <a  (t-Od/]  dx  d.; 
J(j  J()  JQ  10  r  '  n 


9  o  o 

’n  h  Jg  Mun  *  \  *  wn>  dx  d* 


(34) 


1 


9 


For  dn  impact  loading  as  shown  in  Figure  1,  equation  (34)  becomes 
'(•;+. ;  )/a  r  +  .  ft 


it) 


=  '—Z12.  ;-U- 


Lpxix^^)un+p4(x,»,i)vn+pr(x,i,;)wn]sin  .^(t-  ):udxo. 


u.  f  f  oh(u2  ♦  v2  +  w2)  dx  dt 

J0  i0  n  n  n 


(35) 


For  a  concentrated  impact  loading,  equation  (lv4)  becomes 


ijr  *(f+-.  *)/a  r'+* ;  ft 

•u  ...  w.  !  lpx(*.  e,  >)■„  ♦  !>,(*>  5.  >)»„ 


.Ui  ■ 


Hi 


-  (  l)/a  '  Jl! 


+  Pr(x,  $,  x)wrj]  sin  <^(t  -  )  d».  dx  d$ 


a-  ft 


(36) 


m 


!  I  oh(u2  +  v2  +  w2)  dx  d+ 

j(J  j0  n  n  n 


uni  re 


«’x  =  47 


!L/ 

P 


p;  =  4^,, 


"r = 


Solution  for  Impulse 

Consiuer  an  impulse  per  unit  area,  i  (x,  o),  i X(x,  4.)  and  i  (x,  4) 

X  0  I 

acting  on  tne  cylinder  for  an  infinitely  short  time.  The  cylinder  may 
now  be  considered  to  be  vibrating  freely  with  tne  following  initial 
cciiuitions: 

At  t  =  U 


u  =  v  =  w  =  0 


(37) 


11 


i  (x,  ,) 
•u  _  x 


i  ( X,  ;) 


jl-  V(x; 

■it  pTl 


Tin*  displacements  for  free  vii/rations  are  given  as 


u  -  *'  u  (A  cos  w  t  +  b  sin  w  t) 
n;t)  ii  ri  m  n  m 


v  =  v  (;v  cos  u  t  +  B  sin  a.  t) 
‘ ii  n  m  m  m 
m=u 


w  =  V  u  (A  cos  a,  t  +  B,„  sin  u,t) 
m  u  m  m  m 

Substituting  tiie  initial  conditions  (37)  and  (38)  into  equation  (39) 
and  mating  use  of  ortnogonality  yields  the  following 


where 


u  =  u,  B  sin  te  t 
f'1  n  m 
m=U 


v  =  •  sin  t, , t 

n  rn  pl 
m=U 


w  =  V  i<  B  sin  wt 
,  f(I  "  "  m 


fdr  ;  ■. 


_  i  io  I 


(i  u  +  Tv  +  i  w  )  (ix 
'  x  m  0  m  r  nr 


ri  u  t'li  f>,  ,, 

"  I  !  pMu„  +  Vn,  +  WJ  dx  ^ 

Jq  ■[)  w  m  m 


For  a  distributed  impulse,  us  shown  in  Figure.  1,  equation  (41)  becomes 


I 


,)/a  *  +t 
_  *  tn-=  i  )/a  ;  ■.-> 


'  [i  u  *■  i  v  +  i  w  1  dx  d;- 
L  x  m  ^  m  r  mJ 


n 


r2it  f'l  ,  p 

-n  I  ^(u„  +  v‘  +  *£)  dx  di 
n  Ju  Jo 
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(42) 


m  m  m' 


For  a  concentrated  impulse,  equation  (41)  becomes 
m  r('i+t  i)/a  k+ 

-d)  !/  Wa  I. 

.  ‘  (n-t  i)/a  ■»  r- 


lim  r(,,+i  i)/a  cC+t2  _ 

,,  I  Livum  +  i  .vm  +  i  w  ]  dx  d$ 

J (n-t-O/a  Jr-c2  xm  *m  rm 

b  '  2  d) 

fil 

“»  j0  j0  1>h(Un  +  +  V  d*  d* 


(43) 


wnerc 


I 


\  “  4t1"e2 

I. 


t  4t1t2 

I 


1  = 


r  4f- 1  c  ? 

Solutions  for  m  -  q 

For  n  =  0,  equations  (1)  through  (3)  and  (14)  through  (lb)  degenerate 
to  the  following  equations: 


all^  +  paw*  -  Pu"  +  a2px  =  pha2  -7- 

-  a%  - Pv" +  ■  »ha2  If 

-aQx  -  ail  -  pau'  -  Pw"  +  a2pf  =  pha2 


V"  +  “3U'  -!«"+&//!«■ 

-a.,U 1  +  klu’"  -  W,,H-  (M-J-)W}  -  (1  -  u,  )w" 


(44) 

(45) 
(4b) 

(47) 

(48) 


+  (n  -  Pha^ 

+  (Pra  /d)  -  —  w 


(49) 


1 

14 
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Solutions  for  equations  (47)  and  (49)  can  be  determined  from  tu<* 
solutions  given  for  unsyrnmetrical  loading.  For  m  =  0,  equation  (30) 
becomes 


pua^r 

c.  no 


‘3’  "  k’ 


0 


2 

V  0 


a3A  -  kA‘ 


kA^  +  (1  -  c^)x^ 

.  .  .  ,  pha2  2 
k  +  1  *  T  u„o 


=  0 


(50) 


or 


wne  re 


0  4  2 

'  +9^  -  g2*  -  g3  =  o 


Q1  =  ^(1  -  V  +  •'(^jp-  +  2a3)]/k(a5  -  k) 
g2  ■  Cu3  -  a5(l  +  I.)  -  0  -  2a5)(^si)]/k(a5  -  k) 


a3 


=  ohaV  (£haV  .  k  .  ,)/k(  .  k) 


(51) 


Equations  (47)  tn rough  (49)  are  now  uncoupled  and  may  lie  solved 
independently. 

The  solution  for  equation  (48)  is  as  follows: 

For  free  vibrations 


v  =  a,  cos  » ph/U  <*>  ^  +  A,,  sin  ^ph/o  u>  -■ 

I  a  l  a 


(52) 


Trio  ortnogoridlity  conditions  are 


/  2  _  2-. 

("n  v 


pha- 


v  v  dx  = 
n  m 


3V  3V  |  i’ 

y  -■  ■—  -  y  - 11 

11  sx  m  3X  | 


0 


The  forced  vibration  solution  becomes 


U_  [  dx 

n  Jo  n 

T*c  orii«gMlitjr  conditions  fro*  equation  (24)  becone 

{-f  -  ‘£>  ;u  "tV.  *  V.’  ^ 

«u  aw 

L  n*  »i  Ha  a  '  nlww  «x 


♦  %  - p  ir> 

*  '  S1)  *  !l,„  Irf  - 0 


T fee  dynamic  solutions  fron  equations  (33)  and  (34)  are 


(54) 


warn- 


“  *  1-  “on‘x>  V(t) 


n=l 


m  *  l  u  (x)  q  (t) 
ni1  o"  on' 


(55) 


■u“>  ■ 


j„  [px(x-  ')uno  *  pr(x>  ‘Ho*1"  “on(t  *  '>  d'  dx 


I  ^(lli  +  wj"  )  dX 

no  Jq  no  no 


Integral  of  tfie  Square  of  Eigenfunctions 


The  integral  of  tiie  square  of  tiie  eigenfunctions  is  evaluated  from 


')  / 

;  piiLul'ix)  +  v, ,(x)  +  u"(x)]  dx 


=  lin  (-l/dc.  )  7 —  ju  (!«  +  pw  -  P  — — ) 

R  dw  ll'  xm  m  jX 

„  *»„  n 

m  it 

1  i  3v  • i;  .w 

+  v  {;.  -  -  P  — )  -  w  (Q  +  P  — )  +  — !i 

nv  x;m  a  >x  /  nvvxn  3x  xr  -x 

dU  il  ?v  . 

rr  xii  1  n  3x  '  rr  x;n  a  x  ’ 


>w  ^  3w 

+  W  (fj  +  P  — -)  -  il  — — i 

nr  *xn  sx  '  xii  -x  |( 


Illustrative  Example  for  Cylinder  Supported  by  i hi n  Diap'ragn 


For  a  cylinder  su|>ported  by  a  ti.in  diaphragm  the  following  displace¬ 
ments  satisfy  tiie  natural  boundary  conditions  as  derived  from  tnc. 


ortiiogoiial ity  conditions: 


l  l  cos  ri*  cos 

pi  n 


y  ”  V  sin  rie  sin 
u  n  ll" 


\i  =  /.  £  d  cos  me  sin  — - 
1 1,  n 


fo  determine  the  natural  frequencies  and  mode  shapes  tue  determinant 


for  tne  frequency  equation  becomes 


r 


— 

i^V-n  - 
»  < 

1+.  ;;.a, 

~r~  ni - ) 

u 

•3(— )-' 

•I 

I  SI. 

4 

*  t. 

:u 

nii 

•  f 

1 

'  .  1  +  x.  iX\ 

|l  *  (I’ll  ) 

i 

2  /ii-a«2 

-4W  -1(— > 

•  A-O-M - )  1 

♦  <_ 

iV 

j 

1 

i 

,t-lia<~  2 

i 

i 

l 

i 

j 

V  Jnm 

i 

s 

! 

-(21+1-  .  }r~+l 

*♦  t 

1. 

ni 

.  .-ii  d  3 
i  -i-l - ) 

1 

1 

ii 

i 

-  —  /  l 

L  ‘  nii  ! 

(t>8) 


Solutions  for  Forced  Viurations 

Kren  liquations  (33)  and  ( 34 )  tin.’  dynamic  displacements  k.comr 


wliuri 


i 


u 


V 


L 

n 


ri 


:  b 


no 


’  V 


III! 


u 

-  nn 

r. 


cos  m;  (cos  --)qn|J(tj 

sin  n*(uin  ~-)(inn(t) 

cos  ni(sin  ~)q  (t) 
>  nn 


V<1> 


r'-  t 

;  j  j  [p*(x,  >)U  cos  me  cos  ^ 

*o  ht  '0 


+  p<?(x,  >)VW  sin  13$  sin  ^ 


♦  Pr(x,  ♦,  0Unn  cos  m«  sin  J 


sin  ^(t  -  >)  dA  dx  d$ 


“ml  I  ph(Uw  cos2|r^  cos2  ITL 
™jo  Jo 


+  v^sin2n*  sin2  ^ 


+  H2  cos2n$  sin2  £J£-)dx  dj 
mn  * 


Trie  ratio  of  the  mode  shape  coefficients  are 


_  nn  _  -CD  +  Bt 


'  w™.  Au  -  M2 

mn 


-  mn  _  ~/\L  +  BC 

"  "  C”  Ad  - 

mn 


wnere 


rniax2  2  ,  pha  2 


-  _  /rind \C  , 

- - j(_r>  ■"'!* 


1  +  v  _/iiiiak 

i,  =  — - —  m( — - — ) 


r  —  /Una  \  i/l  *  i!/rii'3\  ,  I  /ii’taxS 

c  =  t3(— )  -  kC — 2 — )m  (— )  +  M— ) 


■  i  -  (naa,2  plia^  2 

'  '  "‘4nl  •  Jl(~>  +  iT“  mr 


.  _  ./3  ~  v\„/niTa\2 

t  =  -..^m  -  l;( — ^ — )ni( — r — ) 


For  m  =  () 


"on(lJ 


r2"  r  f‘ 


'U  JO  'U 


[px(x,  V,  A)Uon  COS  -T-+  Pr(x,  <)1*0||  sin  ™ j 


sin  u on ( t  -  a)  d-  dx  d$ 


(,>  dTi  n  j?  (U2  +  VJ-’ 

on  '  on  on 


>\  i 
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For  ra  f  0 


rZi  r9.  ft 

v(ci  =  i  u\(x» *)u, 

'  ; o  -o  Ju  1 


U cos cos  1TL 


♦  l\(x,  +•  ')V  sin  m*  sin  ~ 


♦  Pr(x,  aJW^  cos  m  sin  ~] 


sin  u  (t  -  a)  dA  dx 

■w*  t  <ui + v™,  *  £> 


Solutions  for  Kadial  Impact 


a.  Unit  Step  Loaded  Distributed  over  Finite  Area  (<!c,  U 


4  (t 


-  h-K^] 


*  cos  n<^  cos 


n.x  f  <'  *  cos  “ran1’  i 


l  (‘*Cm(a*  +  6„  + 

v  nm  nm  nm 


„U,fe](™  >H('» 

_ i...  i »  ■  !k  - 

sin  m  sin  —  p  ,(-aT-+  67  yy) 

v  nm'  nm  mn  1 


“  =  TUP 


_ n.x  f  "  -  “5  V‘>  ' 

sin  m  sin  —  -T'(o’r"+e2'  +  1) 

'-iim'nm  pnm 

fell  "  fl]  L.„  n^cl  mie,  cin  Ottx  f1  “  cos  ,unotl 

m  j,  wr  s,n  -f2  s,n — kK0  +  i)i 


^..LU^hTjh^h^h^i 


„  „  (1  -  COS  U)  „t) 

KC0sn*s,n^  <5.v . . ftttt  (w> 

nm'  nm  nm 
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u.  Concentrated  Unit  Step  Load 


c.  Unit  Impulse 

Solutions  for  a  unit  impulse  can  be  found  by  differentiating  with 
respect  to  time  tiie  solutions  for  a  unit  step  function.  A  typical 
displacement  relationship  for  a  concentrated  unit  impulse  is  as  follows: 


;?o 


Triangular  Loading  with  suddenly  nppliod  Valuo  of  Unity,  ami  hctruasi'i:: 

Limarly  u*  /or o  at  Tire,  t 
*  u 


e.  Octangular  Pulse  Loading  with  Suddenly  Appli.u  Value  of  Unity 


and  .juration,  t 


Lxpressions  for  uf  v,  and  w  are  identical  to  those  corresponding  to 
a  triangular  loading  with  tiie  exception  that  Fn(j(t)  and  Fnjn( t)  be 
defined  as  follows : 


Fno(t)  =  (l  -  C0S  “no1* 


Fnn(t)  =  (1  -  cos  “no1*  *  1  c., 


F"oU)  '  ^  H  “„oU  '  V  •  “S  “no1} 

F,«(t)  =  {“*  “>«U  ‘  ‘d*  •  cos  “ran4}  ‘  -  ‘d 
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Equations  of  notion  for  Timoshenko  Theory 


Equations  (  1  )  through  (  3 )  can  i>c  reduced  to  those  presented  by 
Timoshenko  and  Gere  uy  assuming  the  following  conditions: 

a.  The  circumferential  strain  c  .  ew  and  y  are  equal  to  zero 

?  X 

in  calculation  of  X.  and  X  . 

$  Xq> 

b.  i <er.!brdiie  forces  are  not  effected  by  bending  stresses,  nor 
uending  moments  by  imyibrane  stresses. 

Assumptions  (a)  and  (b)  yield  d  *  il  and  IL  =  Ii„..  Assuming 

$X  X$  Qa  X# 

=  (v*  +  w)/a;  u‘,  u"  *  U,  and  v  ~  (u*  +  v‘)/a  =  0,  equations 

V  X?- 

(1)  through  (3)  become: 

+  a,i*x +  pa(v‘* +  w‘) +  =  If  1 


an*  +  ad'  -  a()  -  Pv"  +  a^pA  =  pha^  -r-~r 

\  X$  <J>  0  dt 


C 


-aQ#  -  aQ’  +  all  -  [>a (vT  +  w)  -  PwM  -  a4" =  pha  t t 

•I»  9  r  .iv 


Tlie  rnenbrane  forces  anu  moments  from  equations  (f»)  through  (13)  become 


11.  =  —  (v*  +  u  +  vu') 


il  =  —  (u1  +  vV*  +  vV/) 
x  a 


,1;x^(L^)(u-  tv.) 


i  \r  =  —7-  (“V*  +  w”  +  vw" ) 


il  =  -V  (w"  +  vlf  -  vV ’ ) 

A  a 
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:,,x  =  F-ii  -  W-v) 


‘W  =  F  (1  “  v)(w"  "  v'> 


Substitution  of  equations  (71)  into  equations  (tia)  -  (70)  yields  the 
following: 

u"  *  (  “■■p)vl‘  +  u"  +  vw'  +  q^(v*‘  +  w') 

+  d  P*  -  pha2  u;’u 

+  ~ - D ~W  (72) 

(^-y-^)u'*  +  v“  +  (l^-)vM  +  w*  -  k[w"‘  +  w***] 

+  k[(l  -  v)v"  +  v]  -  q^v"  (73) 

<ju'  +  v*  +  w  +  k[w'*+  2w“"  +  w****]  -  k[v**‘+  (2  -  v)v"'] 

+  q2w"  +  ^  (w  +  w)  -  =  "£jpi  0  (74) 

Equations  (72)  -  (74)  may  be  written 

u"  +  u"  +  ^v'*  +  e^w'  +  a2  (75) 

(~^”~)u 1  ’  +  83v  +  b4v"  +  w*  -  k(w"*  +  w— ) 

+  !^=pi^^ 

L  HFHF  (76) 

\>U'  +  V*  +  ebw  +  k(w"u  +  2w,'•’  +  w’”’)  -  k[v‘“  +  (2  -  v)v"#] 

♦  q2w"  ♦  q,ir  -  tpi  .  -£{jii  0  (77) 


where 


+  (h  = 


2i 


Data  for  Illustrative  Fxar.pl  e  for  Unit  Radial  Impulse 


n  =  1.2  inches 

a  =  60  inches 

i  -  24  incnes 

-,  =  12  inches 

=  2  inches 

«2  *  2  inches 

>  =  ()  radians 

j  =  0.33333  n  =  1  -  30 

n  =  0  -  29 

Pc(F1ugge)  =  S0a3*8b5  PS1  Pc(Tinoshenko)  *  5071*  °3  psi 


■‘t S-  j-+4>*uh 
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ijUL  i 

-HUi  *  .iV.&JUC»i..:ii  I’IIUSU.1  UCi  hrtLQlftiiClLb  1..  bUULl.lU  nUOL 

■  aivrjdil  *  > U.V I 

!*a«eun:  n/2a  -  O.UI  t /2a  =  0.2 

S»j*U  1 1  |1>J  .«Mc:  l!  -  I,  !l  *  9 

!’r-’^sur>_:  V  = 


rlugge's  liieory 

!acli»li»g  filial  Itacriia  .!;:gleccii'»j  axial  Inertia 


p/l*c 

f. 

i 

f, 

4L 

f3 

a  ’ 

fi 

u 

.,17.04 

6UJ/.24 

bl9.7>: 

0..’ 

4o2.*»7 

3929. 7ii 

6003.23 

4b4.91 

3922.72 

b/99.21 

402. b2 

o.«* 

327.03 

i9b.b4 

b79b.  19 

32:;.  74 

O.a 

23 1.2b 

390«^.bb 

6791.  t:s 

232.46 

1.0 
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T/tLLt  11 

I.HLU  UF  HYiMlSlAllL  PiU-SSUKt  Oil  FkLQULMtlLS  IN  BUCKLING  IKJUL 
•  laterial :  steel 

Parameters:  li/2a  =  U.Ol  t/2a  =  0.2 

Luii. ling  iiode:  n  =  1,  i.i  =  9 
Luii.ling  Pressure:  jO/1.793 


Timoshenko's  Theory 

including  /uial  inertia  neglecting  Axial  Inertia 


P/P 

c 

fi 

f, 

C. 

f3 

5 

0 

u  l/.Ulj 

3937.01 

0007.29 

519.02 

0.2 

'(Cl-. 4/ 

3933.90 

0007.90 

404.94 

0.4 

400. 5 i 

3930.95 

0000.52 

402.65 

0.6 

327. 0«i 

3927.92 

0009.14 

320.76 

0 « iJ 

231.24 

3924.90 

0009.75 

232.47 
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T.lbLL  III 

LFFLCi  UF  ISVUHUSTmi  1C  PiitbSUKL  Uil  FKEQUUJCItS  III  bUUXIilb  ;K){iL 
llaceriai :  steel 

Parameters:  ii/2a  =  u.l  i/2a  =  ;* 
buckling  ; iotiu :  n  *  1,  m  =  2 
buckling  Pressure,  =  68811.43 


Flugge's  Theory 


P/P 

c 

Including  axial  Inerti 

fl  f2 

a 

f3 

Neglecting  axial 

fl 

0 

889.93 

6634.11 

12765.31 

997.21 

0.2 

796.00 

bbl3.09 

12750.83 

891.93 

0.4 

689.51 

6593.61 

12736.34 

772.44 

0.6 

5b3.U6 

6573.27 

12721.83 

630.69 

0.8 

398.19 

6552.86 

12707.32 

445.97 

1.0 

— - 
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iaoLI.  IV 

:F»..u  Ul  hYurlioiit'i  lL  PliSSliUL  U.J  FiiOjUlllCILS  III  OUUCLIJiu 

<..Luri.il;  lluuI 

I’.ir  ti  k.iuri:  h/2a  =  0.1  ,:/2a  =  ti 

L.Jti.  iil.cj  I  jjliu:  n  =  i ,  |<|  =  2 
uMt.  nay  Pruusure,  I*  =  00350. <1 

liujsiiunLo's  Theory 

inclu.ii.ig  Midi  Inertia  neglecting  Midi  lnerlid 


iVI*. 

f. 

f , 

t 

i 

L- 

3 

fl 

u 

000  •  l >2 

0023.01 

12011.59 

995.52 

0.2 

793.12 

0023.57 

12010.07 

890.42 

0.4 

ooo.9j 

ou2 3. 34 

12000.54 

771.13 

CJ.li 

500.90 

0023.12 

1200/.03 

629.03 

u.;. 

3%.  71 

0022.09 

12005.51 

445.21 
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TA8LL  V 

CUiiP/iKlSOO  UF  FkLQULWCILS  FOk  VAKlUUS  Him.  SiUvHF  !> 
Flugye’s  Theory 


Parameters:  h/^fa  =  0.01  i/ia  =  0.2 

n  =  1 

P/P^  =  0  P/P .  =  0.0 


1 

fl 

f2 

f3 

f, 

f2 

U 

0 

072. 40 

2079.08 

4471.33 

037.79 

2072.29 

4467.7/ 

1 

060.22 

2001.76 

4003.18 

029.10 

2094.2/ 

4a03.96 

2 

0*1  j.7'j 

2000.71 

4010.37 

000.12 

2008.04 

40 10..  10 

3 

01 9.20 

2770. 30 

4/83.00 

470.08 

2761.32 

4783. bO 

4 

492.44 

290/. Jl 

0021.12 

431.69 

289/. 20 

001 0.00 

5 

471. b3 

3072.0b 

0304.90 

394. b/ 

3061.23 

0298.04 

/ 

O 

40 l.4b 

3201.31 

0032.13 

300.24 

3248.95 

0024.92 

7 

'(•00.20 

34/0.90 

0990.09 

348.42 

3406.0/ 

0987.48 

( • 
iJ 

433.9b 

3t)%.!;7 

6389.04 

347.99 

3680.88 

6379.98 

9 

bl /.04 

3930.33 

0807.24 

300.02 

3919.1!, 

0797.20 

10 

j03. 1  / 

'H08.o7 

7240.89 

400.08 

4109.22 

7234.80 

,  i.  *J*x  WnljMH.  /■ 
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TAliLL  VI 

UftfivUSUlJ  OF  I KLQUkflClLS  FOP.  V/lkIUU5  IM)L  5IIAPLS 


Fiuggo's  Theory 


Paraix-icrs:  h/2a 

=  0.01 

*/2a  =  0.2 

n  =  3 

v/vc  -  o 

P/Pc  =0.5 

u 

T'i 

f2 

f3 

fl 

f2 

f3 

0 

1099. A 

7730.74 

13404.75 

1807.77 

7716.86 

13392.16 

1 

1902.1 1 

7/45.76 

13416.82 

1810.44 

7723.81 

13404.20 

£ 

1910.69 

7/60.78 

13452.99 

1818.46 

7744.65 

13440.26 

I926.U4 

7801.68 

13513.04 

1831.89 

7779.25 

13500.14 

4 

1945.26 

7860.26 

13596.67 

1850.85 

7827.42 

13583.54 

1) 

1971.41 

7912.26 

13703.44 

1875.46 

7888.90 

13690.02 

l. 

2003.66 

7987.35 

13832.83 

1905.85 

7963.35 

13819.04 

7 

2042.11 

80/5.16 

13984.19 

1942.21 

8050.42 

13969.98 

i  ? 

2006.89 

81/5.24 

14156.82 

1984.69 

8149.68 

14142.14 

9 

c 1 30. 1 3 

8287.1b 

14349.96 

2033.47 

8260.67 

14334.76 

10 

2196.92 

8410.40 
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14547.01 
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T/iLLL  VII 

iMVlb&SM  ui  i  iJtQUtnCH.9  Rii;  V/tidUUS  •ii.ii.  SltoPLS 


Flugge's  Tliuoiy 


P.iranexers :  li/2a 

=  U.01 

i/2d 

=  0.2 

ii  =  9 

V/Vc  -  0 

P/P  =0.9 

lil 

T1 

T2 

f3 

fl 

f2 

f3 

0 

bU9i  .(.3 

12097.91 

22339.04 

4990.41 

12061.43 

22310.02 
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b094.0b 
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22347.00 

9001.96 

I2{5b9.99 

22326.04 
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c. 
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12914.67 
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12070.09 

22347.70 
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9120.61 
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22404.99 

9026.04 

12093.00 

22303. 7 9 

*» 
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12964.  <54 
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12927.79 

22434.12 

9 

bl 72.1b 
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9297.97 
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22776.92 
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b  392 .  00 
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10 
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iaLLL  VI i I 

t.FrLii  uh  JiYuHUjI.'.HL  PULSSl  ±  OH  miJWUWIT/J.  lU.QULHCY 
l’.,r.u»:  turj :  li/2d  =  0.01  t/2d  =  0.2  P  =  50i>3.Hbb 


IVPC 

»»• 

n 

f 

0 

1. 

l) 

461.4b 

0.2 

1, 

7 

422.42 

0.4 

1. 

7 

3/4.72 

O.b 

1, 

» 

313.02 

o.<; 

1. 

a 

230.72 

u.% 

1. 

9 

11b.b3 

0.90 

1. 

9 

73.13 

1.00 

1, 

9 

0.00 

***+>■  ••  >.  . 
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T/fcLt  IX 

LKFLCT  OF  iiYUuOSTnTlC  FlitSSUhi  UN  HlUUtP.  HUIQUEHCILS 
Fluggc's  Theory 
Parameters:  ii/2a  =  0.01  i/2  a  =  0.2 


m  =  10 

P/I’c  =  0  P/Pc  =  0.5 


n 

f2 

f3 

fl 

f2 

f3 

1 

663.17 

4100.57 

7245.09 

400.58 

4169.22 

7234.85 

y 

C. 

1196.5b 

6123.37 

10597.00 

1077.53 

6100.93 

10584.18 

3 

2195.92 

0410.40 

14562.79 

2080.67 

8382.92 

14547.01 

4 

3600.27 

10030.29 

10756.16 

3498.39 

10796.81 

10736.09 

5 

5414.29 

13310.62 

23052.99 

5315.33 

13270.66 

23029.97 
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EFFECT  OF  HYDROSTATIC  PRESSURE  u.',  DYNAMIC 
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OYNAMIC  RESPONSE  FOR  UNIT  IMPULSE  WITH  IN-PLANE  INERTIA 
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EFFECT  OF  HYDROSTATIC  PRESSURE  AND  LOADING  AREA  ON  DYNAMIC  RESPONSE  FOR  UNIT  IMPULSE 
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arlous  Values 


Figure  <*.  Longitudinal  Strain  versus  Time  for 


?>1  Inertia 


Circumferential  Strain  versus  Time  for 


Conclusions 


Large  hydrostatic  pressures  and  small  variations  of  impact  area 
greatly  affect  the  dynamic  response  of  deep  submersible  hulls  st Ejected 
to  a  localized  impact  loading. 

For  free  vibrations  deep  hydrostatic  pressures  reduce  the  lower 
frequencies  substantially  while  the  higher  frequencies  are  not  appre¬ 
ciably  affected.  Hydrostatic  pressures  in  the  neighborhood  of  50  percen 
of  the  buckling  pressure  can  reduce  the  fundamental  frequencies  by  30 
percent,  while  the  higher  frequencies,  especially  the  second  and  third 
frequencies  of  the  n,  m  mode  will  have  no  appreciable  change. 

Comparison  of  frequencies  with  the  Flugge  and  Timoshenko  theories 
show  good  agreement  as  illustrated  in  Tables  I  and  II. 

For  forced  vibrations  as  illustrated  by  a  localized  unit  impulse, 
the  following  conclusions  can  be  made: 

a.  Deep  hydrostatic  pressures  have  predominantly  large  effects  on 
longitudinal  displacements  and  strains.  Consequently  the  longi¬ 
tudinal  stresses,  a  ,  will  be  more  sensitive  to  change  while  the 
circumferential  strains  and  stresses  will  increase  moderately. 

b.  Shearing  stresses  experience  moderate  increases  and  are  very  small 
in  magnitude. 

c.  Radial  displacements  and  response  times  will  have  considerable 
increases  as  shown  in  Figure  3. 

d.  Small  changes  ir.  the  area  of  loading  have  tremendous  influence  on 
displacements  and  stresses  as  shown  in  Table  XIV. 
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e.  Comparison  of  theories  indicates  the  following: 

(1)  The  greatest  discrepancy  occurs  in  longitudinal  displacements 
and  strains. 

(2)  Within  the  area  of  impact,  stresses,  radial  and  circumferen¬ 
tial  displacements  have  good  agreement,  while  those  outside 
the  area  of  impact  can  have  large  discrepancies. 

(3)  A  good  estimate  of  stresses,  radial  and  circumferential  dis¬ 
placements  within  the  area  of  impact  can  be  found  by  neglecting 
in-plane  inertias. 
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|  1  2  SPONSORING  Ml  LI  TARV  ACTIVITY 


Department  of  Defense  (Project  THEMIS) 


following  Flugge's  exact  derivation  for  the  buckling  of  cylindrical  shells,  the 
equations  of  motion  for  dynamic  loading  of  cylindrical  shells  subjected  to  hydrostatic 
and  axial  pressure  have  been  formulated. 

The  equations  of  motion  are  applicable  for  long,  short,  or  thick  shells,  and  are 
very  useful  in  calculating  deflections  and  stresses  when  the  impact  loads  are  applied 
to  comparatively  small  regions  of  the  shell.  The  normal  mode  theory  was  utilized  to 
provide  dynamic  solutions  for  the  equations  of  motion. 

Solutions  are  also  provided  for  the  Timoshenko- type  theory,  and  comparisons  are 
made  between  the  two  theories  by  considering  and  neglecting  in-plane  inertia  forces. 

Comparison  of  results  is  exemplified  by  a  numerical  example  which  considers  the 
effect  of  hydrostatic  pressure  on  the  dynamic  response  of  a  shell  simply  supported  by 
a  thin  diaphragm  and  subjected  to  a  localized  unit  radial  impulse. 
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